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Abstract. We investigate effects of disorder on the density of states, the single 
particle response function and optical conductivity in multiband superconductors with 
s± symmetry of the order parameter, where s± — > s++ transition may take place. In 
the vicinity of the transition the superconductive gapless regime is realized. It manifests 
itself in anomalies in the above mentioned properties. As a result, intrinsically phase- 
insensitive experimental methods like ARPES, tunneling and terahertz spectroscopy 
may be used for revealing of information about the underlying order parameter 
symmetry. 
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1. Introduction 

The discovery of iron-based superconductors p] (FeSC) put forward experimental 
and theoretical efforts to understand the reason for rather high critical temperatures 
and symmetry of superconducting order parameters in these compounds. These 
studies yielded a comprehensive experimental description of the electronic Fermi surface 
structure, which includes multiple Fermi surface sheets in a good agreement with 
density functional calculations [2]. The Fermi surface of the moderate doped FeSC 
is given by two small hole pockets around T = (0,0) point and two electron pockets 
around M = (ir, ir) point in the folded zone. This band structure suggests strong 
antiferromagnetic fluctuations, which may be a mechanism for electron pairing. In this 
case, the natural order parameter for most of the FeSC is so called s± state, described 
by the nodeless order parameter with different signs for electron and hole-like pockets. 

This model agrees well with the experimentally found for most of the moderately 
doped Fe-based superconductors nodeless character of the order parameter [31 IU EJ E] . 
However a question, whether the order parameter changes its sign by changeover from 
electron-like to hole-like pockets, is still under discussion. Moreover, relative robustness 
of the superconductors against nonmagnetic impurities led to a suggestion that a more 
conventional two-band order parameter a uniform sign change is realized in these 

systems [7j. 

In our previous paper [8] it was demonstrated that not only superconductors 
with order parameter but also s± may be robust against nonmagnetic impurities. 
Therefore the robustness against nonmagnetic impurities can not be considered as a 
strong argument against s± order parameter. Moreover, it was shown, that there are 
two types of s±-superconductors with respect to disorder [8]. In the first one T c goes 
down as disorder is increased, until it vanishes at a critical value of the scattering rate. 
This behavior is similar to the famous case of Abrikosov-Gor'kov magnetic impurities. It 
is widely discussed in the literature. In the second type of s± superconductors T c tends 
to a finite value as disorder is increased [9j [10]; at the same time the gap functions for 
the electron-like and hole-like Fermi surfaces acquire the same signs, i.e. the transition 
occurs from s± to s ++ . 

In the present paper we discuss how the disorder induced transition s± —> s ++ 
can display itself in single particle properties and optical conductivity. It is shown 
that disorder dependence of these characteristics at the transition point is strongly 
nonmonotonic function of the impurity scattering rate and can be easily seen in tunneling 
spectroscopy, ARPES and optics. Therefore systematic study of disorder effects by 
means formally phase insensitive techniques, mentioned above, may provide information 
about the sign of the underlying order parameter in the clean limit. 

The paper is organized as follows. In the section II we discuss the approximation 
used for the calculations. The section III is devoted to single particle properties. 
We discuss how the transition s± — > s ++ manifests itself in ARPES and tunneling 
spectroscopy. In the section III we calculate the two-particle response function and 
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discuss peculiarities which can be seen in optical conductivity in the vicinity of the 
transition point. The article is concluded with discussions. 

2. The formalism 

For the calculations we employ the standard approach of quasiclassical (£ - integrated) 
Green functions in Nambu and band space [TT] : 



with band quasiclassical Green functions 



ga(w) = -tirN a — (2) 



where the f, denote Pauli matrices in Nambu space and N a is the density of states on 
the Fermi level in the band a = a,b (for the sake of simplicity the two band model is 
considered) . 

The function g a is related to the full Green function 

G a (k,u n ) = Qafo + Uk)f3 + * afl (3) 

by the standard procedure of ^-integration g Q (w) = N a f d£ a (k)G a (k,uj). 

The quasiclassical Green functions are obtained by numerical solution of the 
Eliashberg equations: 

dzK»Jz,u)Re— r ( M=L + iT afi (u)-=SM== 




dzKlp(z,uj)Re _ + iT a p(u) — ? = = ^ , (5) 

where the kernels Kt'n(z,u) have the standard from: 



= J dn 



oo 

tanh ^ + coth ^ 
z + Q — oj — iS 



x 



(6) 



For simplicity we use the same normalized spectral function of electron-boson interaction 
B(Q) for all the channels, which is presented in inset Fig. 2. The maximum of the 
spectra is Q s f = 18 meV [5J. The matrix elements Afg are positive for attractive 
interactions and negative for repulsive ones. The symmetry of the order parameter 
in the clean case is determined solely by the off-diagonal matrix elements. The case 
sign\^ b = sign\f a > corresponds to s ++ superconductivity and sign\^ b = sign\f a < 
to s±. The matrix elements A^ have to be positive and are chosen A^ = |A^|. For 
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further calculations we use the same matrix Xf a = 3, \f b = 0.5, X^ b = —0.2, X ba = —0.1 
for s±-case and Af a = 3, X bb = 0.5, X^ b = 0.2, X ba = 0.1 for s ++ -case. The correspondent 
ratio of the densities of states is N a /Nb = 0.5 (see [8]). 

The second terms in the right side of the Eqs.(jHISJ) reflect scattering on impurities. 
In the general case T a p(oj) can be written in the following form: 

r^M = 7», (7) 

where the 7^ are inter- and intra-band impurity scattering rates in the normal state. 
The dynamical part I {00) = 1 in Born approximation (see Appendix). Beyond the Born 
approximation it reads: 

where C a b(u) is the coherence factor. 

C ab (u) = 1 / _ (9) 

7,2 _ . r,2 _ 12 



Note that in the normal state C a b(u) = and I (00) = 1. 

The dimensionless constant ( is related to the inter-band impurity scattering rate 
in the normal state 7^ as 

il = (io) 

The dependence of 7^ and ( on the scattering potential is shown in the Appendix. 



3. Quasiparticle properties 

3.1. DOS in superconductive state. 

Interband scattering is expected to modify the gap functions and the tunneling density of 
states (DOS) in the superconducting state in a multiband superconductor. In the weak 
coupling regime the impurity effects have been discussed in [TJ] within the Born limit and 
extended in [13] to the strong coupling case. In the following we will calculate the gap 
functions, and the superconducting DOS by solving the nonlinear Eliashberg equations 
in the s± and s ++ superconductors for various values of the interband nonmagnetic 
scattering rate, going beyond the Born approximation. 

Total DOS in superconducting state is given by the following expression 

N( U ) = N a (0) Re 2 " — , (11) 

where we have introduced the complex order parameter: 

A a (a; + i5) = u(f> a (u + i8)/Cj a {u3 + i5) = ReA a (u) + i MA a (w). (12) 

The solution for A a (uj) allows the calculation of the current- voltage characteristic /(V) 
and tunnelling conductance Gns(V) — dI NS /dV in the superconducting state of the 
NIS tunneling junction. 
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In contrast to a single band case, where DOS does not depend on nonmagnetic 
impurities in the multiband case A a (u + id) and DOS are strongly dependent on 
interband impurity scattering. 

Figure 1 shows the calculated gap functions A a (u) for the bands a and b for 
different interband impurity scattering rates. One sees in both s± and s ++ cases strong 
non-monotonic frequency dependence of the gap function with the maximums of the 
absolute values around 250cm _1 for a-band 140cm _1 for 6-band, originated from the 
strong electron-boson coupling. Furthermore, the effects of impurity scattering are 
visible as additional structure at low energies comparable to the interband scattering 
rate r a . The most spectacular effect is the impurity-induced sign change of ReA(u) 
at low energies in fe-band in the s± state. This s± — > s ++ transition was predicted 
in Matsubara representation in our earlier paper (8], where its consequences for non- 
vanishing T c vs disorder behavior were discussed. 




Figure 1. Superconducting gap functions for bands a and b at various interband 
impurity scattering rates in s± and s ++ models. The parameters are C ~ 0.2, 
Ibl = ^laa = lab = ^7ba ~ 0-4ra- They corresponds to scattering strength a = 0.5. 
The relation between a and T a to the scattering potential is given in the Appendix. 



Figure 2 shows comparison of DOS in s± and s ++ states for different magnitudes 
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Figure 2. Total density of states for various impurity scattering rates in s± and s++ 
models at low temperature T <C T c q. The parameters are the same as in Fig. 1. Inset 
shows the electron-boson interaction function B(O). 

of the interband scattering rate T a at low temperatures T <C T c0 . In the clean limit, 
one sees two different excitations gaps for the two bands. In accordance with earlier 
calculations for s ++ superconductors [13], the interband impurity scattering mixes the 
pairs in the two bands, so that the states appear in the a-band at the energy range of the 
6-band gap. These states are gradually filled in with increasing scattering rate. At the 
same time the minimal fe-band gap in the DOS raises due to increased mixing to the a- 
band with strong electron-boson coupling. In the s± superconductor, the modifications 
of low-energy DOS with interband impurity scattering is completely different. Due 
to sign change of ReA(u) in the 6-band, gap less region exists in a range of values of 
scattering parameter T around 35 cm -1 , as clearly seen in the left panel in Fig. 2. Such 
gapless behavior manifests itself in optical properties of s± superconductors, as will be 
demonstrated below. 
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3.2. ARPES and the self-energy 

Angle-resolved photo emission spectroscopy (ARPES) probes the photoemission current 
J(k, co), which can be calculated as: 

J(k, co) = co)\ 2 f(co)A a (k, co). 

a 

Here M(k, co) is the dipole matrix element that depends on the initial and final electronic 
states, incident photon energy and polarization, f(co) is the Fermi distribution function 
and 



A a (k,u) = -^-Tr\ln 



-—Im- 



CO n [CO) 



(13) 



is single particle response function. 

In the weak coupling limit the contribution of the electron-boson interaction to 
self-energy S a0 (k, co) (see first terms in l.h.s. of Eqs. SJ E]) vanishes. It means, that in 
the model with isotropic self-energy I]^ b (u) — > 0, T,^ b (co) — > A a (co). Then the single 
particle spectral function takes the form: 



co 



A a (k, co) = —Im- 

7T 



1 + i E r a p/Jco*-Al(co) 

/3=a,b 



D 



with 



D = ?M + ( E - w2 + A * h - 2 * E r < 

\/3=a,6 / /3=a,f> 



[w 2 + A a (w)A /3 (a;)] 



In the gaped regime A a (k, cj) vanishes below A^ but in the gapless one for 6-band is 
the same as in the normal state: 



co 



Ab(k, co) = —Im 

7T 



1 + i J2 r b/3 /\co\ 

/3=a,b 



2 ' 



V P=a,b 

The quasi-particle spectral function A(,(k.,u) given by Eq. (113ft for 6-band is shown 
in Fig.3. In this case the behavior of A&(k, co) at small co and £ reflects the existence 
of well-defined energy gap. In contrast to that, the function A&(k, to) in the regime 
of s± — > s ++ transition shows no gap, as seen from Fig.4. With further increase of 
scattering rate T a , when s ++ state is realized, in the 6-band energy gap appears again. 
Therefore, ARPES measurements at various impurity concentrations may provide useful 
tool to distinguish underlying pairing symmetry of superconducting state in pnictides. 





100 200 300 400 



&(k) 



Figure 4. The quasi-particle spectral function AbQt,cj) for 6-band with a small gap 
in the gapless regime (T a = 40cm -1 ). The parameters are the same as in Fig. 1. 
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4. Optical conductivity 

The optical conductivity in the London (local, q = 0) limit in a-b plane is given by 

a{u) = ^^JL^uy^riu, (14) 

a 

where H a (oj) is an analytical continuation to the real frequency axis of the polarization 
operator (see, e.g. Refs. [113, [IE], [H], [18], [19]) 



n«(w) = |mT^n a (a;^,i/ m )| 



Wm=HJ+iO+ 



a = a, b is the band index. 



ILvfo;) = / du) 



, f tanh(^) 
\ D R 

tanh (^) — tanh (^) 
~D~ a 



u R Q R + 0^02 



Q R Q R 



tanh (^) 



D J 



1 - 



U) A U)+ + 0^0+ 



AnA 
+ 



Q A Q 



1 - 



UJ A 0J R + (f) A (/) R 

AnR 



Q A Q 



(15) 



where 



and 



D a = VK)M0?) 2 - VR A ) 2 -(0-) 2 , 



w± = u' ± w/2, and the index R(A) corresponds to the retarded (advanced) brunch of 
the complex function F R ( A ' = ReF ± HmF ( the band index a is omitted) . 
In the normal state the conductivity is 



Co' 



pi 



imui 



00 tanh((z + u)/2T) - tanh(^/2T) 



u} a (;z + to) — u a (z) 



If the dominant contribution to the quasiparticle damping comes from the impurity 
scattering, it reduces to the Drude formula: 



Oa(^,T) 



pi 



1 



4tt 7a ^ - iu 



with 77* = 7 ab + 7 aa . 

In the Fig. [5] we demonstrate the impact of disorder on the optical conductivity 
Rea(u). In the clean limit one sees Rea a (u) = for u < 2A a . With increase the 
impurity scattering rate the region Reab{oj) = for the band b decreases and the peak 
above 2A& becomes sharper. It is clearly seen that in the vicinity of transition from 
the s± to s ++ state (r o ~ 35cm -1 ), the conventional Drude- response characteristic 
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in the weak for a normal metal-state is realized. The origin of this effect is gapless 
nature of superconductivity near the impurity-induced s± —> s ++ transition. With 
further increase of the impurity scattering rate, the optical conductivity recovers gapped- 
like behavior, but with smaller gap. It is strikingly different from the behavior of 
superconductors with s ++ order parameter (Fig. [5jb), where values of two gapes tends to 
merge at the limit of infinite impurity scattering rate. This reentrant behavior of optical 
conductivity with concentration of nonmagnetic impurities may serve as unambiguous 
indication for the s± order parameter symmetry. 

Another important characteristic of the superconducting state is the real part of the 
electromagnetic kernel (polarization operator) which is related to the imaginary part of 
optical conductivity Ima(u) (see Eq JT%|) . Fig. 6 shows the frequency dependence of Re 
n(oj) for s± and s ++ models for various interband scattering rates. One can see that 
in the s ++ case dips at cu = 2A a (u) occur for nonzero scattering, in accordance with 
previous calculations done for single-band superconductors [20]. Further, an interesting 
peculiarity is seen in the response of the fe-band in the s± state: the dip position is a 
nonmonotonic function of interband scattering rate and the dip vanishes completely in 
the gapless regime corresponding to the s± —> s ++ transition. 

Magnetic field penetration depth Xl(T) in the local (London) limit in a-b plane is 
related to Imaiui) in zero frequency limit 

1/\\{T) = V lim Antu Im cr a (u, q = 0, T)/c 2 , (16) 

a=a,b 

where c is the velocity of light. If we neglect strong-coupling effects (or, more generally, 
Fermi-liquid effects) then for a clean uniform superconductor at T = we have the 



relation X L = c/uj p i, where uj p ^ a = y 8ire 2 (N a (0)vf a VF a } is the plasma frequency in 
different bands. 

Partial contributions to the magnetic field penetration depth can be written as 

c 4m-o Z a ( u , T) ^ - A|( W , T) [AJ( W , T) - u,*] 

Here the points uj g (T) are determined by the condition for the density of states in the 
band 

ReN(co < u g (T)) = 0. 

For superconductors with gap nodes as well as for T > 0: oj 9 (T) = (see, [2~T]). 

Peculiarities of the penetration depth in the crossover regime from s± to s ++ state 
have been discussed earlier [8]. 



5. Conclusions 



We have studied the effects of the impurity-induced s± —> s ++ transition in the density 
of states, the single particle response function and optical conductivity in a multiband 
superconductors with s± symmetry of the order parameter. It has been shown that 
smaller gap vanishes in the vicinity of this transition, leading to gapless nature of 
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Figure 5. Real part of the optical conductivity Rea(uj) for various values of T a . The 
parameters are the same as in Fig. 1. 
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Figure 6. Real part of the polarization operator for various values of T a - The 
parameters are the same as in Fig. 1. 



photoemission and tunneling spectra. In optical response, the s± — > s ++ transition 
leads to "restoring " of the "Drude"-like frequency dependence of Rea(u). We have also 
found interesting anomalies in the real part of polarization operator, with reentrance 
behavior of the dip-like structure at u — 2A a (c<;) as a function of interband scattering 
rate. This effect leads to non-monotonic behavior in the magnetic field penetration 
depth clS 8b function of the impurity concentration. 

We want to stress that systematic study of the impact of disorder on the single- 
particle response function and optical conductivity may give an information about 
underlying symmetry of the superconductive order parameter. 
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Appendix 



The impurity part of the self energy is obtained by analytic continuation of the 
correspondent self-energy on Matsubara frequencies Ti l ^ p (iuj n ) = (iu n )ro + 

(^ra)^3 derived in [8]. The last is taken in the T - matrix approximation: 

± imp (oj n ) = n tmp ti + tig(ou n )£ mp (iu n ), (A.l) 

where U = U ® % and n imp is the impurity concentration and the scattering potential: 



U 



U aa Ub a 
Uab Ubb 



The quasiclassical Nambu Green's functions on Matsubara frequencies are: 
mN a u an 7rN a (f) an 

90a = , = , 9la = j = ■ ( A -2) 

V &ln + 4>ln V + 4> 2 an 

Solutions of Eqs. (lA.lj) for T^ p ^ and Ea™ p ^ are 

yim P (0)_ _ ^aa-(detU) 2 (^ 2 b -^) » fl „x 
^aa \ 90b + n imp — r gob, l A --Jj 



V<"*(i)-_„ ^ 2 a -(detU) 2 (^-^J _ U ab U ba 

^aa "imp n / i/la "-imp n / \ i/lfej \ ry -^ : ) 

L> \UJ n ) D \UJ n ) 

where 

D(u n ) = 1 - (g 2 0a - g 2 la ) U 2 aa - (g 2 0b - g 2 lb ) U 2 b + (g 2 0a - g 2 la ) {g 2 b - g 2 lb ) (det U) 2 

-2U ab U ba {g 0a gob - g± a g\b) ■ (A.5) 

The analytical continuation on the real axis leads to the following expression: 

and 

0s (cj) 



_^ " up\~y y 

o,6 yj U>}(u) - (ffyiu) 



with r Qi g(a;) given by Eq.([7j). 

The dimensionless constant ( is convenient to express in terms of dimensionless 
scattering potentials u a p = irU a pNp and d — u aa u bb — u ab u ba . Then it has the following 
compact form: 

C = U abUba /. fi N 

^ (d-iy + {u aa + u bb y 1 • ! 
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Normal state impurity scattering rate reads: 

W>imp d -\- U aa 

nN a (d- l) 2 + (u aa + u bb 



N n imp d + U aa m w\ 

7aa = -TTTj 1» , /- , - ( A - 7 ) 



and for a ^ /3 

7q/? = vriy Q (d-l)2 + ( Maa + n fefe )2- 1 • j 

The Born approximation corresponds to <C 1. Then up to quadratic terms in u one 
gets: 



and 



r ab (a;) « 7^ « n^nNbU^. 



It is worth to introduce the parameters cr = u a bUb a / (l + u a bUba) and T a = ni mp ^- = 
n imp^NbU a bUb a {)- — cr). The parameter a is used as an indicator of the strength of the 
impurity scattering. In the Born approximation a — > 0, while in the opposite unitary 
limit a — > 1. 
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